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ABSTRACT 

In view of recent efforts to probe the physical conditions in the pulsar current sheet, we revisit the standard solution 
that describes the main elements of the ideal force-free pulsar magnetosphere. The simple physical requirement 
that the electric current contained in the current layer consists of the local electric charge moving outward at close 
to the speed of light yields a new solution for the pulsar magnetosphere everywhere that is ideal force-free except in 
the current layer. The main elements of the new solution are as follows: (1) the pulsar spindown rate of the aligned 
rotator is 23% larger than that of the orthogonal vacuum rotator; (2) only 60% of the magnetic flux that crosses the 
light cylinder opens up to infinity; (3) the electric current closes along the other 40%, which gradually converges 
to the equator; (4) this transfers 40% of the total pulsar spindown energy flux in the equatorial current sheet, which 
is then dissipated in the acceleration of particles and in high-energy electromagnetic radiation; and (5) there is no 
separatrix current layer. Our solution is a minimum free-parameter solution in that the equatorial current layer is 
electrostatically supported against collapse and thus does not require a thermal particle population. In this respect, 
it is one more step toward the development of a new standard solution. We discuss the implications for intermittent 
pulsars and long-duration gamma-ray bursts. We conclude that the physical conditions in the equatorial current 
layer determine the global structure of the pulsar magnetosphere. 
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1. INTRODUCTION 

The current theory of pulsars is founded on the assumption 
that the magnetosphere of a rotating magnetic neutron star will 
be filled with an electron-positron plasma that everywhere sat- 
isfies the force-free and ideal MHD conditions pE + J x B = 0 
and E ■ B = 0, respectively (Goldreich & Julian 1969). Here, E 
and B are the magnetospheric electric and magnetic fields, and 
p = V ■ E and J are the magnetospheric electric charge and cur- 
rent densities as measured by an inertial (nonrotating) observer. 
Based on this assumption, it has been shown that magnetic field 
lines cannot remain closed beyond the so-called light cylinder 
radius Rlc = c/Q where LI is the neutron star angular veloc- 
ity, and that open magnetic field lines carry the electromagnetic 
energy responsible for the pulsar spindown. Contopoulos et al. 
(1999; hereafter CKF) first showed that the Goldreich-Julian 
requirements cannot be implemented everywhere in the magne- 
tosphere. When one tries to obtain a solution that fills all space 
(and not one with regions of avoidance as in Lovelace et al. 
2006), one finds that it necessarily involves an equatorial cur- 
rent sheet that provides closure for the large-scale electric cur- 
rent that flows from the star to infinity and back. In particular, 
the equatorial current sheet first discovered by CKF continues 
all the way to the stellar surface through two separatrix cur- 
rent sheets that meet at the so-called Y-point where the region 
of closed field lines ends on the light cylinder. This has now 
become the “standard” picture of the pulsar magnetosphere. 

Unfortunately, the CKF solution (and all subsequent improve- 
ments and generalizations of it in ideal MHD) cannot say any- 
thing about how it is implemented microphysically. In other 
words, the solution for E and B requires a uniquely determined 
distribution of p and J, but there is no provision in the solution 
for where the required charges come from and how they are 
channeled around in the magnetosphere. This problem is most 
severe along the equatorial magnetospheric current sheet in the 
interior of which ideal MHD fails and dissipation of electromag- 


netic energy most probably takes place. Therefore, in order to 
understand how the ideal MHD global solution is implemented 
in practice, we need to focus on the physical conditions in the 
equatorial current sheet where the physical requirements are 
most severe. 

Indeed, there has been strong recent interest in pulsar current 
sheet microphysics (Coroniti 1990; Lyubarsky & Kirk 2001; 
Uzdensky et al. 2011; Uzdensky & Spitkovsky 2014; Cerutti 
et al. 2013). All these studies focus on current sheet micro- 
physics independent of the global structure of the magneto- 
sphere. In a series of three articles (Contopoulos 2007a, 2007b, 
2007c), we emphasized the strong coupling between the global 
magnetic field topology and reconnection in the equatorial cur- 
rent sheet. We showed that the nonideal MHD condition in the 
current sheet may dramatically modify the global electric cur- 
rent distribution and consequently the global magnetospheric 
solution as well. These results are confirmed in recent efforts 
to move away from ideal MHD with the introduction of mag- 
netospheric dissipation through various ad hoc prescriptions 
for a finite plasma conductivity (Kalapotharakos et al. 2012; 
Li et al. 2012b). 

In view of the recent efforts to determine the physical con- 
ditions in the pulsar current layer and motivated by our earlier 
2007 work, we believe that there is a need to reconsider the 
main elements of the ideal force-free pulsar magnetosphere. We 
realized that by imposing a simple microphysical requirement 
in the equatorial current sheet we uniquely obtain a new pulsar 
magnetosphere solution that may be used as the new standard 
for calculating pulsar spindown, magnetospheric particle accel- 
eration, and magnetospheric radiation emission. 

2. THE OLD STANDARD SOLUTION 

In the present work we consider only the axisymmetric case. 
We will work in a cylindrical system of coordinates ( R,<j> , z ) 
centered on the pulsar and aligned with its rotation axis with 


1 


The Astrophysical Journal, 781:46 (5pp), 2014 January 20 


Contopoulos, Kalapotharakos, & Kazanas 


4> an ignorable coordinate. We find it convenient to rescale our 
radial coordinate to x = R/Rlc- The magnetic axis may be 
aligned/counteraligned with the rotation axis, in which case 
we will call the pulsar an aligned/counteraligned rotator. We 
also introduce the function 'P = v F(x, z) as the magnetic flux 
contained inside cylindrical radius x at position (x, z). Let us 
begin with a short review of the main microphysical issues of the 
CKF solution where current sheets are treated as dissipationless 
tangential discontinuities. 

1 . The distribution of poloidal electric current 


IcKFm 


("( 2 -*L) iff « iv. 

lo if V F > v l / opcn 


along open magnetic field lines is obtained as an eigenfunc- 
tion of the problem by requiring that the solution crosses the 
light cylinder smoothly. Here, v P opL . n = 1.23 v Pdi po ie, where 
dipole = 7Tr l #*/ R\.c, is the polar value of the neutron 

star magnetic field, and r* is its radius. This contains an 
equatorial current sheet beyond the light cylinder (x > 1) 
with a surface charge density 


ct(x, 0) = 2 E z (x, 0 + ) 

= 2xB r (x, 0 + ). (1) 


cr is positive/negative, therefore, the current sheet consists 
mostly of outflowing positrons/electrons. Their origin may 
be related to pair formation at the Y-point. Note that 


3. Open field lines that cross the zero space charge (null) 
surface contain a distributed outflowing/inflowing electric 
current consisting mostly of inflowing electrons/positrons 
inside the null surface and outflowing positrons/electrons 
outside. Their origin may be related to pair cascades in the 
so-called outer gaps. 

4. A small amount of open field lines cross the null surface 
twice and contain a distributed inflowing/outflowing elec- 
tric current. The origin of its associated charge carriers is 
obscure. 

5. Most open field lines around the pole do not cross the null 
surface. Those contain a distributed inflowing/outflowing 
electric current consisting mostly of outflowing electrons/ 
positrons. Their origin may be related to pair cascades on 
the pulsar polar cap or in the outer gaps. 

The most important part of the magnetosphere is the equato- 
rial current sheet that guarantees the global poloidal electric cur- 
rent closure. The current sheet experiences a pinching force due 
to its reversing surface magnetic field components 5^(x, 0 + ) = 
—B<p{x, 0~) and B r (x 1 0~) = — B R (x, 0~). The new element 
that we want to emphasize here is that the current sheet is also 
charged because the electric field component perpendicular to it 
also reverses across it, i.e., E z (x, 0 + ) = — E z (x, 0~). Therefore, 
the vertical electrostatic self-repulsion of the current sheet sur- 
face charge density a contributes to its vertical support against 
the magnetic pinching force. Note that, in contrast to the solar 
corona, there is no “guide field” in the pulsar reversing equato- 
rial current sheet as some recent investigations assume (Cerutti 
et al. 2013). 


er(l + , 0) = 0 (2) 

right outside the magnetospheric Y-point on the light 
cylinder, since B R ( 1 + , 0) = 0 at that location. The fact 
that the equatorial current sheet does not only carry an 
electric current but is also charged has not been adequately 
emphasized in previous studies. As we will see below, 
this has some very important consequences for the global 
structure of the pulsar magnetosphere. 

2. The equatorial current sheet closes onto the star through 
two current sheets along the separatrix between open and 
closed field lines. The separatrix current sheets are charged 
with a surface charge density 

er(x, z) = E(x, z + ) - E{x, z~) 

= — xB p (x, z + )- 


3. THE EQUATORIAL CURRENT SHEET 

Let us now study the small-scale structure of the current sheet 
and assume that it has a thickness 2h <£ R. In the current sheet 


interior defined as —h ^ z ^ h, in general 

B R {x,z)=jB R {x,h), (4) 

h 

B z (x, z) = B z (x, h), (5) 

B,p(x, z) — jBpix, h), (6) 

n 

E r (x, z) = E r (x, h) = —xB z (x, h ), (7) 


1 + 


[2/( v t , opeD )/(Rc)] 2 
X 2 [ 1 — x 2 ]Z? 2 (x, Z + ) 



(3) 


E,(x, z) = f E z (x, h) = xB r (x, z), (8) 

h 

E< j> {x,z) = 0. (9) 


cr is negative/positive, therefore, the separatrix current 
sheets consist mostly of inflowing electrons/positrons. 4 
Their origin too may be related to pair formation at the 
Y-point (Contopoulos 2009). Here, B p is the poloidal 
magnetic field. 


Note that Lr(x,0) = /L(x, 0) = 0 in CKF. In general, 
E r (x, 0) = xB z (x, 0) > B(x , 0). We then have 


Jcs 


Bc/,(x, h)c~ 

K + 

h 


Br(x, h)c 
h 


0 , 


( 10 ) 


4 For the history of this discovery, we must acknowledge that although CKF 
did indeed discover the magnetospheric current sheet associated with the 
poloidal magnetospheric current, it was the referee of their article, Professor 
Leo Mestel, who pointed out the presence of tangential discontinuities at the 
two separatrices. 


(T cs = 2 E z (x, h) = 2 xBr(x, h) , (11) 


Pcs 


E z (x, h) 
h 


xBr(x, h ) 
h 


( 12 ) 
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Figure 1. (a) Distribution of electric current 47r/New/<^ / dipoie as a f unc ti° n of 
'F/'Fdipoie along open field lines. T* 0 pen = l-4'F<iipoie- fNewC^open) = 0. The 
solid/dashed/dotted lines indicate new/CKF/monopole solution, respectively. 
The dash-dotted line marks the value of T'max i = 0.8 5 dipole = O.bOT'open 
beyond which the return magnetospheric current flows, (b) Distribution of 
Poynting flux normalized to £ckf- The solid/dashed lines indicate new/CKF 
solution, respectively. The dash-dotted line marks the value of T'max i beyond 
which the magnetospheric Poynting flux gradually enters the current sheet. 


The vertical force experienced by the current sheet plasma is 
equal to 


F z (x , z)cs — pF z + 


JrB,/, — 


[(x 2 - 1 )B 2 r (x, h) - B 2 (x, h)]z 


h 2 


03 ) 


In the old standard pulsar magnetosphere solutions, zF, is 
negative, i.e., the inward magnetic pinching force dominates. 
Most importantly though, at the points where the equatorial 
current sheet meets the light cylinder, 


F Z (\,Z) = 


Bl(hh)z 

h 2 


(14) 


In standard solutions, B ^(\ , h) = 2/(T , ope n)/(Rc) ^ 0 and, 
therefore, F z ( I , z) ^ 0. An equivalent way to say that the 
magnetic pinching force in the current sheet dominates is to 
realize that the electric current in the current sheet in the old 
standard solution is spacelike, i.e., that 


(B 2 (x,h) + Bl(x,h)?' 2 c 
h 

E,(x, h)c 

= Pcsf. 

h 


05 ) 


This picture requires extra charge carriers of the opposite sign 
(electrons/positrons) in excess of the ones needed to account 
for p cs . It also requires a “hot” population of particles to support 
the current sheet against its own magnetic pinching. 

A very simple and natural simplification to our problem is to 
assume that the electric current in the current sheet is everywhere 


null, or at least only weakly spacelike. In that case, Equation (13) 
yields F z ~ 0. In other words, a null equatorial current sheet 
is electrostatically supported by its own charge and requires no 
additional “hot” population of particles. The above requirement 
amounts to a new equatorial boundary condition given implicitly 
through 


2Q IQV(x,h)) 


— x(x 2 — l) 1/2 ZJ s (x, h) 


( 16 ) 


or equivalently 


Q.IC¥(x,h)) _ tJ1 lU/2 

o — {X J ) 

7 TC Z 


aT(x,z) 


dz 


Z=h 


(17) 


instead of the “standard” CKF equatorial condition v P(x, h) = 
F(x, 0) = ¥(1,0) = ¥ open . Note that 2/(¥(x, h))/(Rc) = 
xB^ix, h). 


4. THE NEW STANDARD SOLUTION 


As described in CKF and Timokhin (2006), we solve the 
pulsar equation numerically by iteratively relaxing both the 
distribution of ¥(x, z) everywhere in the half space z > 0 (given 
the boundary conditions along the axis, infinity, and the equator), 
and the distribution of magnetospheric electric current / Ob) 
that allows for a smooth and continuous crossing of the light 
cylinder. The main difference with CKF is the new equatorial 
boundary condition (Equation (17)) that is implemented during 
the iteration to evolve the equatorial distribution ¥(x , 0) for the 
next iteration step. 

Our results are shown in Figures 1-3. In Figure 1(a) we show 
the distribution of poloidal electric current that we obtained 
along field lines that cross the light cylinder (we will still call 
them “open”). This can be fitted with the analytic expression 





T open / \ T open / 


0.4 


if M' < d' opcn 

if^^^open ’ 


Here, T , ope n = 1 .4 v |-' ( | lpll | e . For comparison, we also show 
the CKF current distribution and its corresponding monopole 
distribution. In Figure 1(b) we show the distribution of Poynting 
flux in the pulsar magnetosphere, starting from the axis (*P = 0) 
up to the last open field line (*P = v P opcn ). For comparison, we 
also show the CKF Poynting flux distribution. 

The main elements of the new solution are as follows. (1) The 
pulsar spindown rate is equal to 


E = 0.82 


Bl&rl 

4c 3 


( 18 ) 


i.e., 82% of the CKF spindown or 23% higher than that of the 
orthogonal vacuum rotator. (2) Sixty percent of the magnetic 
flux that crosses the light cylinder, M' lnax i, carries the one part of 
the magnetospheric current that consists of electrons/positrons 
that outflow to infinity. (3) The electric current closes along the 
other 40% that consists of positrons/electrons that gradually 
converge to the equator. (4) This transfers 40% of the total 
pulsar spindown energy flux in the equatorial current sheet, 
which is then dissipated in the acceleration of particles and 
in high-energy electromagnetic radiation within a few times 
the light cylinder radius. 5 (5) There is no separatrix current 


5 Ninety percent of the dissipated energy is dissipated within six light 
cylinder radii. 
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Figure 2. New structure of the pulsar magnetosphere with a null equatorial 
current sheet. Shown are magnetic flux (*F) surfaces plotted in intervals of 
0.1 dipole- The dashed line indicates T'maxi = O-SS'Fdipde (corresponds to 
the boundary between outgoing and ingoing magnetospheric current, and 
asymptotically approaches the null line — i.e., the boundary between positive 
and negative space charge — and the equator as R — > oo). The dash-dotted line 
indicates 4 / 0 pen = 1 •4'Fdipoie- The thick gray line indicates the equatorial current 
sheet. 

sheet. (6) Pressure balance across the equatorial current layer is 
guaranteed by the electrostatic repulsion of the charge contained 
in it. Our solution is a minimum free-parameter solution in that 
the equatorial current sheet does not require a thermal particle 
population for vertical support against collapse. 

The most important element of that solution is (5), namely 
the fact that 

W^open) = 0 (19) 

(Equation (17) for x = 1 ), and B,j, ( \ , h) = 0. In other words, the 
electric current closes within the field lines that cross the light 
cylinder (we will still call them “open”), and thus there is no 
need for a current sheet along the separatrix between open and 
closed field lines. The equatorial current sheet is established to 
provide closure for the global pulsar electric circuit. As we said, 
only magnetic field lines that carry the return current consisting 
of outflowing positrons/electrons enter the current sheet and 
then outflow along the current sheet to infinity. The main 
magnetospheric current consists of electrons/positrons that 
outflow along the other 60% of magnetic field lines that originate 
closer to the polar cap axis. It is interesting that Equation (19) 
naturally resolves a complication in the CKF solution, namely 
that B p (x, 0)] = 27 C KF( l I / open)^/c 2 [x 2 - l) 1 / 2 !^^ 

oo (Uzdensky 2003; Contopoulos 2009). 



Figure 3. Same as Figure 2 but on a larger scale. 

Note the qualitative similarity between the solution shown 
in Figures 2 and 3 and the solution obtained in Figure 5 of 
Contopoulos (2007c), which was based on analogous qualitative 
global arguments. In that solution, however, there was no 
provision for the ultimate opening up of field lines that carry 
the main magnetospheric current, and therefore the associated 
spindown rates obtained in that work are different. 

5. DISCUSSION 

We have presented a solution that we hope will become the 
new standard for calculating pulsar spindown, magnetospheric 
particle acceleration, and magnetospheric radiation emission. It 
is the first global solution that self-consistently takes care of the 
equatorial current sheet microphysical conditions in the limit of 
a null current. In that limit, the current sheet is self-supported 
electrostatically, and in that respect it is a “minimum parame- 
ter” solution, as it does not require an extra population of hot 
positrons/electrons. We must acknowledge that the electromag- 
netic (Poynting) energy that is deposited in the current sheet 
is dissipated in the acceleration of positrons/electrons by the 
Er = xB- > B / 0 electric field component in the current 
sheet (Equation (7)). Obviously, some fraction of those par- 
ticles is expected to form a hot population that will partially 
contribute to the electrostatic support of the current sheet. In 
that respect, our present solution must be considered as one 
more step toward the development of a new standard solution. 
However, it is interesting that since particle acceleration pro- 
ceeds beyond the light cylinder, any hot population will 
also form beyond the light cylinder. The absence of a hot 
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population at the Y-point justifies one of the main results of the 
present work, namely the elimination of the separatrix current 
sheet. 

The electrons and positrons that populate the pulsar magneto- 
sphere may naturally be produced through pair production at the 
crossings of the zero space charge (null) surface, the so-called 
outer gaps. The electrons/positrons that move inward continue 
to flow through the polar cap and subsequently outflow in the 
main part of the polar cap, whereas the positrons/electrons that 
move outward reach the equatorial current sheet beyond the light 
cylinder (these are the ones that absorb the electromagnetic en- 
ergy dissipated in the current sheet). Our minimum parameter 
solution is physically appealing in that it guarantees that pair 
production is not needed anywhere else in the magnetosphere. 

Our solution addresses a fundamental issue that other ideal 
force-free solutions cannot address, namely, what is the fraction 
of the released electromagnetic energy that is deposited in the 
current sheet and dissipated in the acceleration of current sheet 
positrons/electrons. Our numerical integration yields 40%, 6 and 
it is expected that some significant fraction of it will be seen 
as high-energy pulsar radiation. We have thus found the main 
region of magnetospheric dissipation: the equatorial current 
sheet within a few light cylinder radii beyond the light cylinder. 

Our results may also find application in the phenomenon of 
intermittent pulsars that are observed to spin down faster during 
their ON state than during their OFF state. If we associate the 
ON state with our present solution, the spindown rate of the ON 
state is only 82% of that of the CKF solution. We can naively 
generalize Equation (18) in three-dimensional in analogy to 
Spitkovsky (2006) as 

B 2 Q. 4 r 6 

£on(«) ~ 0.82 * * (1 + sin 2 a) , (20) 

4c ^ 

where a is the inclination angle between the magnetic pole 
and the rotation axis. If we also associate the OFF state with 
the Li et al. (2012a) solution with nearly vacuum conditions 
on the open field lines and nearly ideal force-free conditions 
on the closed field lines (where plasma remains trapped in the 
absence of pair production), then we obtain a ratio of respective 
spindown rates £on/£off 1-2.5 for a likely range of pulsar 
inclination angles of 90°-30°. This may naturally account for 
observed values between 1.5 and 2.5. 


Our solution may also find application in the study of long- 
duration gamma-ray bursts (GRBs) where the central energy 
source is believed to be a newly formed maximally rotating 
stellar mass black hole that loses energy through electrodynamic 
processes very similar to the ones that take place in the 
pulsar magnetosphere (Blandford & Znajek 1977; Komissarov 
& Barkov 2009; Contopoulos et al. 2014). In particular, a 
maximally rotating 10 M G black hole can lose up to about 10 55 
erg of rotational kinetic energy in the form of Poynting flux 
(Contopoulos et al. 2014). However, it is still not clear what 
fraction of that spindown energy will be radiated in the form of 
high-energy (X-ray, y-ray) radiation. Our new standard solution 
yields a significant fraction of the spindown luminosity. If that 
result is confirmed, it will allow us to use long-duration GRBs 
as standard candles in Cosmology. 

We conclude that the physical conditions in the equatorial 
current layer significantly modify the global structure of the 
pulsar magnetosphere. 

This work was supported by the General Secretariat for 
Research and Technology of Greece and the European Social 
Fund in the framework of Action “Excellence.” 
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